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Abstract 

General first- and higher-order intertwining relations between non-stationary 
one-dimensional Schrodinger operators are introduced. For the first-order case it is 
shown that the intertwining relations imply some hidden symmetry which in turn 
results in a i?-separation of variables. The Fokker-Planck and diffusion equation are 
briefly considered. Second-order intertwining operators are also discussed within a 
general approach. However, due to its complicated structure only particular solu- 
tions are given in some detail. 

1 Introduction 

It is a basic fact that many physical phenomena are mathematically described by solutions 
of linear ordinary and partial differential equations. For example, the dynamics of a 
classical system may be characterized by Newton's equation or equivalently by the Euler- 
Lagrange equation. Another example is the dynamics of a non-relativistic quantum system 
which is governed by a linear differential equation, the well-known Schrodinger equation. 
Hence, exact solutions of such kind of linear differential equations are of basic interest 
and, therefore, much effort has been made during the last centuries to find solutions of 
problems being of the form 

= o , (l) 

where L denotes some linear differential operator. Here ip is the wanted solution of ([!]) 
in some function space associated with given initial and/or boundary conditions. Many 
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methods have been developed to find such solutions. This paper is concerned with one of 
them based on an assumed intertwining relation 

C x l = 1C 2 (2) 

between two physically relevant differential operators Cx, £2 and the so-called intertwining 
operator 1, which is also assumed to act linearly. This intertwining relation allows to 
construct a solution for, say, Cx if a solution of £2 is known. To be more precise, let us 
assume that ip2 is a solution of £2^2 = 0, then due to the intertwining relation ipi = lip 2 is 
a solution of Cxipx = 0. Note that in addition we have to assume that 1^2 does not belong to 
the kernel of X, i.e. Iip2 ^ 0. Multiplying @ from left and right with the inverse operator 
X -1 we have l~ x C\ = /^X -1 and, therefore, one can also start with a given solution ipi 
and obtain a new solution ip2 = l~ x ^x ^ 0. However, in many cases the intertwining 
operator is a differential operator. Therefore, its inverse is a rather complicated integral 
operator and of less use. This situation, however, changes if in addition we assume that 
both operators £i j2 are self-adjoint on some common function space, C\ 2 — £1,2- As a 
consequence one has the adjoint intertwining relation 

l ] C l = C 2 1 ] , (3) 

which in addition provides us with the symmetry operators XX"'" and X^X obeying |1J 

[£ 1; XXt] = 0, [£ 2 ,XtX]=0. (4) 

The above connection between solutions of two differential equations have, to our 
knowledge, been applied first by Darboux Pj to Sturm-Liouville problems. That is, £1,2 
belong to a class of second-order differential operators and may, for example, represent 
the Hamiltonian of a quantum mechanical degree of freedom in one dimension. This fact 
may be the reason why most applications of such intertwining relations (sometimes also 
called Darboux transformations) have been made for eigenvalue problems of stationary 
Schrodinger Hamiltonians in one dimension. Here, Cx,2 represent a pair of Schrodinger 
Hamiltonians intertwined by a first-order time-independent differential operator. This 
pair of Hamiltonians together with the intertwining operator (also called supercharge) 
form the basis of Witten's model of supersymmetric (SUSY) quantum mechanics ||. It 
is also closely related to the factorization method of Schrodinger M and Infeld and Hull 
||, and still finds application in the construction of new potentials giving rise to exactly 
solvable Schrodinger-eigenvalue problems. See, for example, || [7|]. This approach has 
also been generalized || |9| to higher-dimensional systems where matrix Hamiltonians 
necessarily participate in the intertwining relations. Very recently, intertwining relations 
with supercharges of second (and higher) order in derivatives were investigated for one- 



dimensional [jlO| , lilf and two-dimensional systems [[I] as well as for matrix Hamiltonians 
in one dimension [JT2]j . In the one- dimensional case it was shown that irreducible second- 
order transformations, these are those which cannot be represented by two consecutive 
standard (first-order) Darboux transformations, do exist. In two dimensions second-order 
irreducible transformations allow to avoid intermediate matrix Hamiltonians. Both, for 
one-dimensional matrix and two-dimensional scalar quantum systems non-trivial symme- 
try operators @ were constructed Jj], and led to integrability of the corresponding 
dynamical systems. 

In this paper we will consider intertwining relations of the non- stationary Schrodinger 
operator 

S[V] =idt + 8£-V(x,t) . (5) 
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Here d t = d/dt and d x = d/dx denote the partial derivatives with respect to time and 
position. Throughout this paper we will denote these derivatives, if applied to some 
function /, by a dot and prime, respectively. Hence, we use the notation f(x,t) = 
(d t f)(x,t) and f'(x,t) = (d x f)(x,t). It is also evident that we are going to use units such 
that Planck's constant 2nh and the mass m are given by h = 2m = 1. 

Replacing the general operator L by the Schrodinger operator fl5|) the corresponding 
equation (|I|) reads S[V]ip(x, t) = and represents the non- stationary Schrodinger equation 
for a one-dimensional quantum system under the influence of an external time-dependent 
potential V : R x R — > R, (x, t) i— > V(x, t). A first investigation of this type of intertwining 
relations is due to Bagrov and Samsonov |L3| who considered a time-dependent first-order 
intertwining operator. See also the review |T^] where in addition a reducible second-order 
intertwining of non- stationary Schrodinger operators is considered. 

In the next section we will reconsider the first-order intertwining and show that, due 
to the i?-separation of variables, it provides no basic new results besides those already 
known from the stationary case of time-independent potentials. This observation has 
also recently been anticipated by Finkel et al. |15|]. However, we additionally show that 



this fact can be related to some underlying symmetry structure. We will also briefly 
discuss the time-dependent Fokker-Planck equation in this context and find that this 
case does not lead to a trivial result in the above sense. In Section 3 we will present 
particular solutions for stationary potentials intertwined by a non-stationary, i.e. time- 



dependent, second-order operator. We show that the reducible cases (as studied in ||T4|| ) 
are only very special cases of this class. In Section 4 we finally discuss the intertwining 
of non- stationary Schrodinger operators by non-stationary second-order operators. As 
special cases we consider the intertwining of the Schrodinger operators corresponding to a 
non- stationary Hamiltonian with that for a stationary one as well as with that of a time- 
dependent harmonic oscillator. In the concluding remarks in Section 5 we discuss some 
aspects omitted in the main text such as zero modes, domain question and generalizations 
to complex potentials. 



2 Intertwining by first-order operators 

In this section we will consider the most general first-order intertwining relations for the 
Schrodinger operator as defined in (|5]) and the Fokker-Planck operator -F[i7], 

F[U] = -d t + d 2 x + d x U'(x, t) = -d t + d 2 x + U'{x, t)d x + U"(x, t) . (6) 

In the latter case, the relation (^) becomes the well-known Fokker-Planck equation 
characterizing the diffusion of a macroscopic (one-dimensional) degree of freedom in an 
external time- dependent drift potential U : R x R — > R, (x, t) i— ► U(x, t) with a diffusion 
constant set equal to unity. For both cases we start with the most general intertwining 
operator of first order given by 

q + {x,t) =Z (x,t)d t + £ 1 (x,t)d x + £ 2 (x,t) (7) 



with, in general, complex-valued functions £o;£i an d £2- Note that in contrast to fl3 
also allow for a first-order operator in dt- 



we 



3 



2.1 The Schrodinger operator and separation of variables 

For the above defined Schrodinger operator (|5]) the most general first-order intertwining 
relation reads 

SIV^ = q + S[V 2 ] , (8) 

where the functions & (i = 0, 1, 2) and the potentials V\$ are clearly not independent 
of each other. In fact, we will now consider the question: Which £'s lead to real-valued 
potentials and how are the latter connected with the them? 

Inserting the explicit forms of the Schrodinger operators (|5]) and the intertwining 
operator (0) into relation (p]) and equating the coefficients of identical (partial) differential 
operators we immediately find that £ an d £i may only depend on time, i.e. £q = = £[. 
The assumption that £ does not vanish identically then leads to the conclusion that 
the potential difference V\ — V 2 does depend only on time. This is a rather trivial caseQ 
and, therefore, we will set £o = without loss of generality.^] Making now the following 
reparametrizations = e 1 ^ p(t) and £ 2 (x,t) = e 1 ^' p(t)x'{x, t) with (3 : R — > R, 

p : R -> M+ and x ■ R x R -> C we find 

V x {x, t) = X ' 2 {x, t) + X "(x, t) - i X (x, t) + a{t) - Pit) + ip{t)/p{t) , 

(91 

V 2 {x, t) = x'Hx, t) - X "{x, t) - i X {x, t) + a(t) , 

where a is some time-dependent complex- valued integration constant. Following the argu- 
mentation made in footnote [1] we may set (3 = without loss of generality. Furthermore, 
we may also set a = because it can always be absorbed in x y i a f ne substitution 
X ^ X ~ i J dta. Hence, we are left with 

V 1 {x J t)=x' 2 {x,t) + x"{x,t)-i X {x,t) + ip{t)/p{t) , 

V 2 (x, t) = X ,2 (x, t) - X "(x, t) - ix(x, t). ' ' ' 

Here the so-called superpotential x is s f m n °f arbitrary as the potentials are assumed to 
be real-valued. This can, for example, be achieved by assuming a stationary real-valued 
superpotential. This, however, leads to the standard stationary supersymmetric quantum 
mechanics discussed to a large extend during the last two decades 0. Therefore, we will 
consider a complex-valued superpotential 

X(x,t) = h(x,t) +ig(x,t) (11) 

with real- valued functions h and g ^ 0. The reality condition Im V\ = Im V 2 = then 
leads to 

2g" + p/p = 0, 2tig'-g"-h = 0, (12) 
which can easily be integrated to 

9{x,t) = ~^&x 2 + \p{t)p{t)x + 1 {t) , 

x 4 PW 2 (13) 

h(x,t) = -\np(t) + K(x/p(t)+fJL(tf) , 



^ote that from V 2 (x,t) - V x {x,t) = f(t) and S[Vx]ipi = follows that S[V 2 ]ip2 = S[Vi + f]tp 2 = 
has the solution ip2{x,t) — exp{z J Q dr /(r)}^i(a;,t). That is, knowing a normalized solution of the 
time-dependent Schrodinger equation for V\ one immediately has a solution for V 2 , too, which is also 
normalized. Complex- valued /'s are not allowed as both, V\ and V 2 , are assumed to be real- valued. 

2 A posteriori, this justifies the ansatz made in fbl. 
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where p and 7 are arbitrary real-valued function of time and K is an arbitrary real- valued 
function of the variable y = x/p + p. In terms of these functions the final form of the two 
partner potentials is 

(14) 

and the intertwining operator reads 

q + (x, t) = p{t)8 x + K'(x/p(t) + //(*)) - ^ (p(t)x - p 2 (t)p(t)) . (15) 
It is obvious that with appropriate choices for K, p, p and 7 at least one of the two 



potentials (14) can be made stationary. We will not investigate this aspect here but 
rather like to discuss a more general property of these quantum systems, the so-called 
i?-separation of variables. 

In doing so, let us demonstrate that the non-stationary Schrodinger equation 

S[V ia }^ 2 {x, t) = (id t + d 2 x - V h2 )^ 2 (x, t) = (16) 



with potentials as given in ([14]) , which is equivalent to the intertwining (HD, admits a 
separation of variables. In fact, under the transformation 

y = x/p(t) + p(t) , 1)2 (y, t) = v^e^Vyfx, *) (17) 

this Schrodinger equation becomes [T7| 

ip 2 (t)d t <f> 1>2 (y, t) = [-d 2 y + K' 2 {y) ± K"{y)] <f> li2 (y, t) , (18) 

which is obviously separable in y and t. Hence, the solutions of the original problem 



(16|) have the general form ip(x,t) = R(y,t)Y(y)T(t) which is known as the ^-separation 
of variables [|T^]. In other words, for any pair of Schrodinger operators ^[Vi^], which 
admits a first-order intertwining relation ([8]) there exists a transformation ([H]) to some 
new coordinate in which the potentials become stationary and form the well-known pair of 
stationary SUSY partner potentials 0. As a consequence nothing essentially new emerges 
from the first-order intertwining of the non-stationary Schrodinger operator. This fact has 
also recently been observed by Finkel et al. [13]. Note that the transformation associated 



with the special case p(t) = 1 and p{t) = vt with constant velocity v corresponds to the 
Galilei transformation. See, for example, the textbook by Galindo and Pascual JIT| . 

Here let us add that this i?-separation of variables can be related to the existence of 
symmetry operators for the two Schrodinger operators in questions. First, we note that 
with the results above one can directly verify the adjoint intertwining relation 

q-S^} = S[V 2 }q- (19) 

where 

q- = ( g +)t = _p( t ) dx + K'(x/p(t) + pit)) + % - (p(t)x - p 2 (t)p(t)) . (20) 

In other words, in analogy to the discussion of @ the Schrodinger operator S[V] and 
in particular its time-derivative idt can be viewed as self-adjoint operators on a suitable 
domain, (id t Y = id t . Then from (|]) and (|19|) immediately follows 

S[V 1 ],q + q-]=Q, \s[V 2 },q-q + ] =0, (21) 
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which is a special case of (^). This discussion suggests that the reverse conclusion may 
also be true. In fact we are able to prove the following 

Theorem: Any Schrodinger operator S[V] (with real- valued potential V) having a self- 
adjoint symmetry operator R, obeying [,S[V],-R] = and being quadratic in d x , admits 
the i?-separation of variables for the corresponding time-dependent Schrodinger equation 
and an intertwining relation of first order with its superpartner. 

For the proof of this theorem we start with the most general quadratic and self-adjoint 
ansatz for the symmetry operator, 

R(x, t) = -u(t)c% + i{8(x, t), d x } + ((x, t) , (22) 

obeying the symmetry condition 

[S[V], R] = [idt + dl- V{x, t), R(x, t)] = 0. (23) 

Here, 5 : R x R — > R and (:MxR->l are arbitrary real functions of position and time, 
whereas uo : R — > R + is, without loss of generality, assumed to be an arbitrary strictly 
positive function of time only. The coefficient functions of the symmetry operator ( |2"2"| ) 
are clearly related with the potential V due to the symmetry condition (|23"D giving rise to 

6(x i t) = \u(t)x + lu{t) i C(M) = $(.) + -f^* 2 + 4^*, (24) 
where <3> is an arbitrary real function of the variable 



dT Zmr) (25) 



and v is another arbitrary real function of time. In terms of these functions the potential 
must take the form 

i . 1 /..,. U 2 (t)\ n 1 /.,, V>(t)u(t)\ $(z) , . 

V(x,t) = TT[U}(t) Tt)x 7X\v(t) - ,\ ]X+ ■ 26 

We note that the i?-separation of variables for the Schrodinger equation S[V]ip(x,t) = 
becomes explicit by making the change of variables 

t) = ^ ( t) exp (H * + J*| , - \[ *■ ) } ^ () _ (27) 
which results in 

(iuj{t)d t + d\ + $(z)) (j)(z, t) = 0. (28) 

Let us remark that the symmetry operator R of the form ( p2|) with coefficient functions 
related by (^1) is factorizable, i.e. R = q + q~ > with q~ = (q + Y . The supercharge q + 
has again the form ( I5|) if the following substitutions are made: 

p(t) - ^(t) , y^z, K' 2 + K"^^. (29) 

Thus we conclude that if some non-stationary Schrodinger operator S[V] has a symme- 
try characterized by a second-order operator of the form (|22]) it admits the /^-separation 
of variables and the superpartner potential of V can also be constructed. The non- 
stationary Schrodinger operator associated with this superpartner has again a positive 
symmetry operator, which in this case is given by q~q + . 
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2.2 The Fokker-Planck equation 

In this subsection we will now consider the most general first-order intertwining relations 
allowing to construct from known solutions of one Fokker-Planck equation solutions of its 
corresponding superpartner. In order to keep as close as possible to the previous discussion 
we will first transform the Fokker-Planck equation into a diffusion (or imaginary-time 
Schrodinger) equation. Indeed, it is well-known |16| that the Fokker-Planck equation 
F[U)P(x,t) = transforms under the substitution P(x,t) = exp{—U(x,t)/2}ip(x,t) into 
the diffusion equation 

D[V)ip(x, t) = {-d t + d 2 x - V{x, t)) if>(x, t) = , (30) 

where the potential V is given in terms of the drift potential U in the following way 

V(x, t) = i U' 2 (x, t) - i U"(x, t)- l ~ U(x, t) . (31) 

Note that the diffusion operator D[V] can be obtained from the corresponding Schrodinger 
operator S[V] via a Wick rotation idt — > —dt- There are, however, essential differences 
from the physical point of view as the solution P of the Fokker-Planck equation should 
represent a probability distribution. As a consequence if)(x,t) = exp{U(x,t)/2}P(x,t) is 
also required to be real and positive. 

We are now interested in finding the most general first-order intertwining operator q + 
as given in (0) obeying the relation 

D\Vi]q + = q + D[V 2 ] <^> F[U X \ e^Ve"^ 2 = e Ul ' 2 q + e~ u * /2 F[U 2 ] • (32) 

Here we closely follow the derivations of the previous subsection and the result is identical 
to that given in ([J) after Wick rotation. Again we may set a = as it can be absorbed 
in x- Furthermore, as noted above, the intertwining operator q + = e l ^p(d x + x') should 
map a positive solution of -D[V2]^2 = into a positive solution of I?[Vi]^i = 0, i.e. 
ipi = q + ip2 > 0. This leads us to the conclusion that (3 must vanish identically and x nas 
to be real-valued. Hence, we are left with 

V 1 {x,t) = x ,2 {x,t) + x"{x,t) + X {x,t)-p{t)/p{t) , 

V 2 (x, t) = X ' 2 (x, t) - X "(x, t) + X (x, t), 

giving rise, cf. (|3~TD , to the following form of the drift potentials 

U 1 (x,t) = 21np(t)-2 X (x,t) , 

(34) 

U 2 (x,t) = 2 X (x,-t) . 

In contrast to our result in the previous subsection we find here that x has to be 
a real-valued function which, however, is a function of the two independent variables 
x and t. Therefore, the Fokker-Planck case can, in contrast to the Schrodinger case, 
not be transformed into a stationary problem. This seems to be related to the fact 
that the diffusion operator is, in contrast to the Schrodinger operator, not self-adjoint. 
Consequently, the adjoint intertwining relation (|3|) and the resulting symmetries @ do 
not exist. Physically, this is due to the time-reversal symmetry, which exists in the case 



of the Schrodinger equation jLJl but not for the Fokker-Planck and diffusion equation. 
This is also explicated in the above result ([54]) which shows that the two SUSY-partner 
drift potentials are essentially related via a time reversal. The additional difference in an 
overall sign between them is already known from the stationary case B. 
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3 Second-order intertwining for stationary potentials: 
particular solutions and symmetry operators 

Recently higher-order generalizations of the SUSY-intertwining relations for stationary 
Schrodinger operators have been investigated to some extend pi, ITT]. In the 1-dimensional 



case it was shown [11] that in general the intertwining of a pair of stationary Schrodinger 
operators by a second-order differential operator cannot be reduced to two consecutive 
transformations with some intermediate self-adjoint Hamiltonian. In the 2-dimensional 
case second-order intertwining operators also allow to intertwine pairs of standard station- 
ary Schrodinger Hamiltonians with scalar potential All of these self-adjoint Hamilto- 
nians have a hidden symmetry, cf. characterized by a differential operator. 

In the remaining part of this paper we will investigate the intertwining of a pair of 
Schrodinger operators S[Vi] and S[V 2 } by second-order intertwining operators of the form: 

q + (x, t) = G(x, t)d 2 x - 2F(x, t)d x + B(x, t) . (35) 

As in the first-order case it can be shown that the inclusion of an additional term being of 
first order in d t leads to the trivial situations where the difference V\ — V 2 depends on the 
time t only, cf. footnote |l|. Furthermore, from the intertwining relation @ with above q + 
one can conclude that the function G may not depend on x and similar to the discussion 
in the previous section it is even possible to exclude a phase. In other words, without loss 
of generality we have G(x, t) = g(t), which should not be confused with g used previously, 
and consider from now on an intertwining operator of the form 

q+(x, t) = g(t)d 2 x - 2F(x, t)d x + B(x, t) . (36) 

We are unable to find the general analytic solutions of the intertwining relation (|8|) 
with q + as given above. Therefore, in this and the next section we will construct only 
some particular solutions of interest. More precisely, in this section we shall look for the 
solutions of the intertwining relation (|8|) for the case where both potentials V\ and V 2 are 
stationary, i.e. do not depend on t. It is evident that one class of solutions is already 
known from ||11|| . Assuming a supercharge q + whose coefficient functions are real and 



do not depend on t, it follows that the corresponding solutions of (J8|) will coincide with 
those of the stationary intertwining relations (— + Vi(x))q + (x) = q + {x){—d^. + V 2 (x)) 
which can be found in JTTJ. Here we are interested in more general solutions of @ with 



operators q + depending manifestly on t. That is, we will search for particular solutions 
of the intertwining relation 

(idt - #i)g + (x,i) = q + (x,md t - H 2 ) (37) 

with standard stationary Hamiltonians H 12 = —d^+Vi^x) but an explicit time-dependent 
intertwining operator. 

3.1 Systems with symmetry operators of third order 

A first suitable ansatz with simple t-dependence in (|36"D is 

q + (x, t) = M + (x) + A(t)a + (x) , (38) 
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where we have set 

M + (x) = d 2 x - 2f{x)d x + b(x) , a + (x) = d x + W(x) . (39) 

Here all functions besides A are considered to be real valuedP]. With this ansatz the 
intertwining relation fl3~T|) results in the condition 

iAa + = H X M + - M + H 2 + A (H x a + - a + H 2 ) . (40) 



Equating the coefficient of identical powers in d x on both sides of (|40f ) one obtains the 
following system of equations: 

iA = 2m + 2mA , 

H 1 M + — M + H 2 = 2ma + , (41) 
H\a + — a + H 2 = 2ma + , 

with real constants m and m. For our further discussion we will consider the two cases 
m^O and m = separately. 

a) Solutions for m ^ 0: 

In this case the first equation in fl41~D immediately leads to A(t) = m§e~' 2%mt — rh/m, 
itiq G K, and results in an intertwining operator of the form 

q + (x, t) = d 2 -( 2 fix) + -)d x + b(x) -—W(x) + m e- 2imt a + (x) . (42) 
V mJ m 

It is obvious that without loss of generality we may set fh = as a non- vanishing rh may 
always be absorbed via a proper redefinition of / and b. As a consequence, the second 
relation in (|4"1~D leads to a second-order intertwining between Hi and H 2 . This has already 



been considered in [11 and it was found that the potentials Vi, V 2 and the function b can 



be expressed in terms of / and two arbitrary real constants a and d: 

f"{x) p(x) d 



Vt, 2 (x) = T 2f(x) + f(x) + 



2f(x) 4P(x) AP(x) 



"(x) f' 2 (x) d 



(43) 



b(x) = -f'{x) + f\x) - i-±-L + ±-±-L + 



2f(x) 4P(x) 4/ 2 (x) 



Finally we note that the last equation in ( flip coincides with the usual first-order inter- 
twining relation of SUSY quantum mechanics between the Hamiltonians Hi and H 2 + 2m. 
Therefore, the potentials can be expressed in terms of the SUSY potential W as follows 
I: 

V x {x) = W 2 (x) +W\x) , V 2 (x) = W 2 (x) - W\x) - 2m . (44) 

Comparing this with (ff^) we conclude that W(x) = —2f(x) — mx and / must satisfy the 
so-called Painleve IV equation: 

f 2 d 
f" = ±— + 6/ 3 + 8mxf 2 + 2{m 2 x 2 -m + a)f + — . (45) 
2 J £J 

As an example we mention here that for d = —a 2 it can be shown that any / obeying 
the generalized Riccati equation /' = — 2f 2 — 2mxf — a is also a solution of (f45|). The 

3 We also assume A ^ 0, as the case A = has already been studied in |Q|. 
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general solution of this generalized Riccati equation and the corresponding potentials 
are discussed in detail in || , which also shows their connection to the most general SUSY 
partners of the harmonic-oscillator potential. Note that cases with d < are reducible 

ones irm . 



We have already discussed in the previous section that the intertwining relation (|37| ) 
and its adjoint give rise to the symmetry operators q + q~ and q~q + for (id t — Hi) and 
(idt — H 2 ), respectively. For the present case of stationary potentials one can use results 
of |lTJ to show that these symmetry operators (defined up to a square polynomial in #1,2) 
have the explicit form 

Rt(x,t) = e 2imt M + (x)a-(x) + e - 2imt a + (x)M- (x) , 

(46) 

R 2 (x,t) = e 2imt a~(x)M + (x) + e~ 2imt M+ (x)a- (x) , 

where a" = (a + Y and M~ = (M + )^. Let us note that these results have also been found 
recently by Fushchych and Nikitin p0[ using a different approach. 



b) Solutions for m = 0: 

In this case it follows from (|4"1|) that A(t) = —2imt where, without loss of generality, we 
have set the integration constant A(0) to zero. The remaining two conditions of (|4~1~D lead 
to 

Vi 2 {x) = W 2 {x) ± W'{x) , f(x) = n-- W(x) , 

1 2 (47) 

b(x) = - (W'(x) - W 2 {x)) - 2nW{x) - mx , 

where W is a solution of the Painleve II equation 

W" = 2W 3 + AmxW + k (48) 

and n, k are real constants. Here we note that an additional integration constant appearing 
on the right-hand side of the last relation in ( fj?D has been set to zero as it can always be 
absorbed via a proper redefinition of the independent variable x. 

Again let us mention that for k = —4rh a particular solution of this equation reads 
W(x) = 1/x giving rise to inverse-squared potentials Vi t z- O n the other hand, for k = 2m 
a solution of the ordinary Riccati equation W = W 2 + k x is also a solution of (fKJ). Note 
that such solutions can be expressed in terms of Bessel functions pl |. 



As in the previous case the intertwining relation (|37|) leads to symmetry operators 



which can also be reduced to third-order differential operators: 

R x = M+M- - H 2 + 2imt(M + a- - a + M~) + Afn 2 t 2 Hi , 

(49) 

R 2 M \T - H 2 + 2imt(a~M + - M~a + ) + Am 2 t 2 H 2 . 

However, in addition to that, because of m = 0, we can even construct another pair of 
third-order operators commuting with (id t — Hip). Indeed, the last two equations of fl4i"D 
can be rearranged as follows 

[Hi, M + a~] = 2mHi , [H 2 , a'M + ] = 2mH 2 (50) 

and directly lead to 

[id t ~ H u iM + a~ + 2mtH x \ = [id t - H 2 , ia~M + + 2mtH 2 ] = , (51) 
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Hence, besides those given in fl49|), we have found an additional hermitean pair of sym- 
metry operators: 

R x = % (M+ar - a + M~) + AmtH 1 = RJ2m , 

(52) 

R 2 = i (a~M + - M-a + ) + AmtH 2 = R 2 /2m . 

Whereas these operators have already been given in [^(J the previous pair fl49|) is a new, 
independent one, which does not commute with (|52|). 



3.2 Systems with symmetry operators of fourth order 



As a second ansatz for the intertwining relation (p7|) we consider an operator with two 
coefficient functions depending manifestly on time: 

q + (x, t) = 9(t)M + (x) + i\(t) x a + (x) , (53) 

where 9 and A are assumed to be real-valued functions and the operators M + and a + are 
of the same form as in (|39l) . For this ansatz the intertwining relation (|37|) leads to 

9 = -2A , \ = P0, 

(54) 

2M+ = xa + H 2 - H x xa + , /3ia + = M + # 2 - H 1 M + . 

The first two equations, containing the real constant (3, can easily be integrated and yield, 
for j3 > 0, oscillatory solutions. Furthermore, with the above conditions the potentials 
Vi t 2 and the functions W and b can exclusively be expressed in terms of the function /: 

W(x) = -2f(x) + - , b(x) = f'(x) + 2f 2 (x) - V 2 (x) + (P/4)x 2 + a , 



x 



2(1 -2c) fix) 13 9 A(c-l)f(x) 
Vx, 2 (x) = T 2f'(x) + Af(x) + '-^ + U 2 + ~ (55) 

Ob O Ob 



x 



dz[f(z)-2zf 2 (z)} +a 



where clq, c and Xq are some real constants and the function / must satisfy the non-linear 
third-order differential equation: 

- 2fV 2 > + Abf + /"' + 4/' 2 + Aff" + 2(5xf -p c + P/2 = (56) 

This equation is certainly too complicate to solve for arbitrary values of the constants (3 
and c. However, for c = and an arbitrary d one can show that all /'s satisfying the 
Riccati equation 

f(x) = 2f 2 (x)-((3/A)x 2 + d (57) 



are also solutions of (pq) . Such solutions correspond to Hamiltonians obeying the inter- 
twining relation H\a + = a + H 2 in addition to (|54]). 

It is straightforward to show that there are no other values^ of the constant c for 
which the Hamiltonians Hi 2 are intertwined by operators of first or second order. Since 



solutions of eq. (|57| ) do not exhaust the solutions of (p6|) we believe that, at least in some 
cases, operators q + of the form fl5"3"D intertwine the Schrodinger operators S^Vi^] but not 
the Hamiltonians H\ 2 . 



1 We exclude here the case = which corresponds to the intertwining of H 1 2 by M + . 
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3.3 Reducible and irreducible second order intertwining. 



The natural question we will now consider concerns the reducibility of the second-order 
intertwining operator ([!6|) to a pair of consecutive first-order operators af and af. We 
again assume that the real potentials V\ y2 = Vi, 2 {x) are stationary but the intermediate 
real potential V may depend on time as well and to which we may always add some 
arbitrary time-dependent function A = A(t). In other words, we consider the following 
chain of Darboux-type transformations between the Schrodinger operators: 

S[V] — > S[V + A] X S[V 2 ) . (58) 

As we have shown in Section 2 a first-order intertwining for S[Vi] and S\V 2 ] implies the 
existence of second-order symmetry operators, that is, 

[S\y 1 ],a+a^] = [S[V 2 ),a^a+) = 0. (59) 

Let us first consider that case where afa± is a trivial symmetry operator, i.e. afa± = 
Hi + const. Then the intertwining operators af being of the form (0) with superpotential 
(PI) can be simplified to 

p(t) = l, g(x,t) = g(t), h{x,t) = h{x). (60) 



Thus the operators af do not depend on time, af = ±d x + h'(x). The same argumentation 

l 2 



applies to the operators af and H 2 . Evidently, in this case the Hamiltonians Hi 2 belong 



to the class considered in |TTJ , where it has been shown that reducible as well as irreducible 
second-order intertwining operators exist. 

In the case of a non-trivial symmetry operator afa~[ the potential V\ is of the form 



see [^0 and references therein) 



Vi{x) = k + hx + k 2 x 2 + 3 —. (61) 

(x + fc 4 ) 

It is evident that this class of potentials does not exhaust all solutions of intertwining 
relations fl37|) with ansatze (|38f) or (|53D for q ± . 

The conclusion drawn from this discussion is, that the intertwining relations inves- 
tigated in this section have both reducible and irreducible solutions of second order in 



4 Second-order intertwining for non-stationary po- 
tentials 

In order to consider the intertwining of Schrodinger operators with a non-stationary po- 
tential, it is useful to simplify the general form (|36|) of the intertwining operator q + a 
little further. Indeed, by inspecting the general intertwining relation (|8|) with q + as given 
in (|36|) one finds the condition 

\mF{x,t)= l -g{t)x + gi{t) (62) 
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with an arbitrary real function g\. After multiplying both sides of (|j) with 



a-W(t) exp \i + ^ - f dt> ^ 

9 [ }e P [ \8g(t) + g(t) J dt g\t>) 

and some algebra we obtain a new intertwining relation, 



id t + 2i 



( gx I 9i 



\±9(t) g 



+ ^\d x -H 1 \q + = q + [id t + 2i(?- + ^\d x -H.. 



4# <? 



(63) 



(64) 



which, of course, is equivalent to @. Here the "new supercharge" and "new Hamiltonians" 
are given by 



q + = gd 2 x -2(ReF)d x + B + -g- 
= gd 2 x -2(Re F)d x + B 



1 f gx 



g\± 



+ 9i 



2iF f gx 



9 V4 



+ 9i 



and 



H\o = Hi o + 



1 ('9 9 2 \ 2 t 1 f 9i9 2<3 



9 2<r 



XT + 



2 \g 



9 



X 



(65) 



(66) 



respectively. 

Now, in a second step we make a non-local transformation [[T7| of the independent 
variables (x, t) 1— > (y, r) with 



A ' 



y = g^l\t)x-2 / dt g^g-*^) 



(67) 



which brings the intertwining relation fl64|) into the form 



z<9 T + ^ - ^ ^ - 2#~2 (Re F)^ + B 



d 2 - 2g-2 (Re F)d y + B ) (id T + d 2 - gV 2 



where 



1 ( '9 9 2 \ 2 , 1 ( 9i9 29 



(6* 



(69) 



,9 2^ 

This shows that, without loss of generality, we could have chosen from the very beginning 
the operator q + to be of the form 



q + (x, t)=d 2 x - 2f(x, t)d x + b(x, t) + ic(x, t) 
with all coefficient functions being real-valued. 



(70) 



4.1 Intertwining of a non- stationary Hamiltonian with a sta- 
tionary one 

In Section 2 we have only briefly mentioned that a first-order intertwining between 
Schodinger operators for a non-stationary potential V\ and a stationary one V 2 is pos- 
sible. Here, because of the absence of the .R-separation of variables, we will reconsider 
this aspect for the case of a second-order intertwining operator of the general form (|70|) 
in more details. 
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The intertwining relation (H) with ansatz (j70|) leads to the following conditions for the 
coefficient functions: 

f = c>, b + c" + Acf' = 0, f»- b >- V > + 4ff' = 0, ( 
c + 2/^'-6"-46/'-V r 2 " = , ^i = y 2 -4f. 1 J 

It may easily be verified that a particular set of solutions of (|7lD can be expressed in 
terms of two arbitrary real- valued functions fi and /o, which depend only on x and t, 
respectively: 



2 h(x)+f (t) " ' y 2 h(x)+f (t) 

v 2 (x) = + ^ (/^/r^) - ^ r - 1 /i 2 (^)) , (72) 

f) = V 2 {x) - Af'(x, t) , b{x, t) = f'(x, t) + 2f(x, t) - V 2 (x) . 

Recall that b, c, f and V\ are functions of x and t whereas V 2 is assumed to be stationary, 
i.e. independent of time t. This latter assumption actually leads to the explicit form 
fo(t) = crexp (Xot) + 5exp (— X t) with arbitrary constants a, 5 and A . For this case, the 
symmetry operators Ri associated with the Schrodinger operators S[Vi\ read 



R x = q+q- = Hi + \ Ag - 2c - Aic'd x - lid' 
R 2 = q-q+ = H 2 2 + \\ 2 , 



(73) 



with Hi = —&l + Vi. Note that Ri explicitly depends on time, whereas R 2 is stationary, 
a fact, which helps to show that the transformed function ipi = q + ip 2 of a normalized 
solution of S\y 2 ]ip 2 (x,t) = is a normalizable solution of S[Vi]ipi(x,t) = 0. Indeed, 
because of the linearity of the Schrodinger operator and the stationarity of V 2 , we may 
set ip 2 (x,t) = ip 2 (x,t; E) = e~ lEt ipE(x) with <£e being an eigenstate of H 2 associated with 
eigenvalue E. Therefore, we have 

fr(x,t;E),Mz,t;E)) = (q + Mx,t;E),<l + M^t;E)) = {E 2 + \ 2 /4)6z E . (74) 

Finally, we mention that in the particular case of a constant potential V 2 (without 
loss of generality we may set V 2 = 0) the function fi has to be a solution of the linear 
fourth-order differential equation with constant coefficients (see eq. ([72])) 

f[" - Xlh = , (75) 

whose general solution may easily be constructed by standard methods and will be an arbi- 
trary linear combination of the fundamental solutions sin(v^Ao x), cos(a/Aox), sinh(\/Ao x ) 
and cosh(\/Ao x). Note that the partner potential Vx(x, t) = —2 [d 2 \og(fx(x) + fo(t))] has 
the same scattering properties as V 2 = 0. That is, it is a reflectionless potential. As 
particular example we may choose fi(x) = cosh(x) and fo(t) = cosh(t) which leads to 
Vi(x,t) = — 2cosh(x) sinh(t)/(cosh(x) + cosh(t)) 2 . 



4.2 Intertwining of a non- stationary Hamiltonian with a time- 
dependent harmonic-oscillator Hamiltonian 

In the remaining part of this section we will now consider the second-order intertwining 
of Schrodinger operators associated with two non-stationary Hamiltonians. In doing so 
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we will limit ourselves to one of the simplest non- stationary Hamiltonians, namely, that 
of a harmonic oscillator with a time-dependent frequency and search for the class of 
non- stationary potentials V\ which are intertwined with this one. 
As a starting point we choose the potential to be of the form 

i (g{t) a 2 {t) , 2 



^-^m-m ]x ' (76) 



which corresponds to the trivial potential V 2 = in as we also set g\ = in the 
following discussion. As a consequence we obtain from ( |68~D the intertwining relation: 

id T + d 2 - U,) (d 2 y - 2f(y, r)d y + B(y, r)) = (d 2 - 2f(y, r)d y + B(y, r)) (z3 T + d 2 ) , 

(77) 

where 

^i(y,r)=^)Vi(y,r), f(y,r) = g- 1 / 2 (t)ReF(x,t) (78) 

are given in terms of the new variables y = x/^fg and r = j^dt' g~~ l {t'). Taking into 
account the relations (^) we have 

^(y,r) = -4(9 y /)(y,r) , S(y,r) = r) + 2/~ 2 (y, r) + ic(y, r) 

/(, x) = 1 c(, r) = l - (79) 

2/i(y) + /o(r) ' ^' 2/ l(y ) + / ( r )' 

where Jo( r ) = 5-e A ° T + 5e~ Af,T and f\ is an arbitrary solution of ( fT5f) with Ao and x 
replaced by Ao and respectively. Returning back to original variables (x, t) and setting 
g(t) = p 2 (t), we find that the original potentials read 



vdx.t) = -M. x 2 - 2 [ai log (/, {p- l {t) x ) + /„ (/ 'dt>- 2 («'; 



*0) 



which may be compared with the result of Bluman and Shtelen |T7|| who derived new 
potentials being related to the free quantum system (in the y-r coordinates) via non-local 
transformations. Here we note that the special case pit) = cos(a;t) leads to the ordinary 
(i.e. time-independent) harmonic oscillator case and the corresponding transformation has 
been called Jackiw transformation.^ 

The present result is also closely related to the special case K l2 — K" = const, discussed 
in Section 2.1 where we have shown its close connection with the R-separation of variables. 
Thus solutions for Vy in (p0|) can be related to a free quantum system characterized by 



5 Concluding remarks 

In this paper we have considered the most general time- dependent first- and second-order 
intertwinings of non- stationary Schrodinger operators. Whereas, in the first-order case a 

5 See, for example, p^, |2j| and p3] where more general transformations are given. 
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complete discussion was possible, in the second-order case only particular examples have 
been discussed. These, in turn, are not as simple as those of the first-order case. 

In the main text we have concentrated on the implications of the intertwining rela- 
tions and did not take into account secondary aspects such as zero modes and domain 
questions, which we are going to discuss now. As already mentioned in the introduction, 
the intertwining may provide a new non-trivial solution if>i = q + ip2 of a given problem 
only in the case when ip2 does not belong to the kernel of the intertwining operator q + . 
Elements of the kernel of q + lead to trivial solutions. On the other hand, S'fVi]^ = 
may have non-trivial solutions which cannot be obtained from solutions of iSf^]^ = 0. 
In this case, ipi necessarily belongs to the kernel of the adjoint intertwining operator q~, 
q~ipi = 0. Therefore, in order to find a complete set of solutions one has to consider also 
the kernel of the intertwining operator and that of the adjoint one. Again the first-order 
case allows a general discussion. The result of our discussion in that case (Section 2) has 
led us to the general form q + = p(d x + h! + ig f ), where h and g are given in ([13]). As q + is 
a first-order differential equation the dimension of its kernel is at most one- dimensional: 

ip2(x, t) oc exp{— h(x, t) — ig(x, t)} . (81) 

The restriction of solutions to the linear space of square-integrable functions then leads 
to the normalization condition 

cfa# 2 (x,t)| 2 = f dyexp{-2K(y)} . (82) 



Hence, the question of zero modes of q + or its adjoint is equivalent to a discussion of broken 
versus unbroken SUSY of the corresponding stationary problem in the ^/-coordinate || . 
The situation in the case of second-order intertwining is a bit more complicated. However, 
it is clear that the corresponding kernels are at most two-dimensional. In the particular 
case studied in Section 4.1 the kernel of q + is obviously empty for a real non- vanishing 
parameter A as in this case the symmetry operator R 2 = (q + Yq + , cf. eq. (|73"|), is strictly 
positive. 

Another aspect which we did not discuss in the main text is concerned with domain 
questions. Indeed, one usually is interested only in square-integrable solutions of the 
Schrodinger equation. There is no guarantee that the intertwining operator maps a square- 
integrable solution into another square-integrable one. This typically happens when the 
superpotential or, more generally, the coefficient functions of the intertwining operator 
become singular M . Hence, the rather general and abstract results presented in the main 
text should be understood in the sense that whenever there appear singular coefficient 
functions a careful analysis of such domain questions is mandatory. For example, the 
case of inverse-square potentials briefly mentioned in Section 3.1 b) is only well-defined if 
the domain (Hilbert space) is defined by the space of square-integrable functions on the 
positive half-line with proper (e.g. Dirichlet) boundary conditions at x = 0. 

Finally, let us mention that it is also possible to extend our approach to complex- 
valued potentials Vi t 2, which recently have attracted some attention ||25| . The case of 
stationary complex potentials with first-order intertwining has already been studied in |7], 
2"E|| . However, because of the absence of the reality conditions, cf. eq. (0), the first-order 



intertwining of such potentials will not lead to the i?-separation of variables. Furthermore, 
the corresponding Schrodinger operators are no longer self-adjoint and, therefore, no 
symmetry operators of the form q + q~ and q~q + exist. This situation is similar to that 
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of the Fokker-Planck and diffusion equation discussed in Section 2.2. For the higher- 
order intertwining the notion of irreducible transformations is lost in the case of complex 
potentials. 
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